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Abstract 

■ We compute nonequilibrium dynamics for classical-statistical SU(2) pure 

q . gauge theory on a lattice. We consider anisotropic initial conditions with 

\ high occupation numbers in the transverse plane on a characteristic scale 

" ~ Qs- This is used to investigate the very early stages of the thermal- 

| . ization process in the context of heavy-ion collisions. We find Weibel or 

"primary" instabilities with growth rates similar to those obtained from pre- 
vious treatments employing anisotropic distributions of hard modes (par- 
ticles) in the weak coupling limit. We observe "secondary" growth rates 
'<-]'■ for higher-momentum modes reaching substantially larger values and we 

Q-f analyse them in terms of resummed loop diagrams beyond the hard-loop ap- 

O proximation. We find that a coarse grained pressure isotropizes "bottom-up" 

with a characteristic inverse rate of y~' - 1 - 2fm/c for coarse graining mo- 
mentum scales of p < 1 GeV choosing an initial energy density for RHIC of 
e - 30GeV/fm 3 . The nonequilibrium spatial Wilson loop is found to exhibit 
. an area law and to become isotropic on a similar time scale. 

J£j '■ 1 Introduction 

(N 



o 



The theoretical understanding of the thermalization process in collision exper- 
iments of heavy nuclei at RHIC (BNL), the LHC (CERN) or the planned FAIR 
facility at GSI is crucial for a successful outcome of the heavy-ion program. Avail- 
able data from RHIC reveals remarkable properties such as rapid apparent ther- 
malization with robust collective phenomena. There seems to be a good descrip- 
tion of a wide range of data using hydrodynamical models, provided one assumes 
a fast approach to local thermal equilibrium in less than about lfm/c [0Q. The 
understanding of these extremely short time scales from QCD provides a chal- 
lenge for theory. It has been pointed out that complete local thermal equilibrium 
may not be necessary for the application of hydrodynamics, and a prethermaliza- 
tion of the equation of state could be sufficient Filamentation or so-called 
Weibel instabilities in anisotropic plasmas have been identified as the paramet- 
rically fastest processes which may help to explain a rapid isotropization of the 
equation of state on time scales much shorter than characteristic thermal equili- 
bration times 
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During the past years several groups investigated plasma instabilities in non- 
equilibrium gauge theories within different frameworks. Extensive studies have 
been performed within the "hard-loop approximation", which is equivalent to the 
collisionless Vlasov theory 0, or solving the effective Wong-Yang-Mills equa- 
tions [J5j . Another important development concerns the solution of classical Yang- 
Mills fields in the presence of static colored source terms using the framework of 
the color glass condensate effective theory in expanding geometries Also 
transport or kinetic equations in the spirit of earlier "bottom-up thermalization" 
scenarios 0, i.e. assuming that no plasma instabilities occur, were performed 
with remarkably short thermalization times O. 

Typically, these studies rely on a sufficiently large separation of scales between 
suitably defined soft and hard momentum scales. This is applicable only at small 
characteristic running gauge coupling and limits many results to the weak cou- 
pling regime. One would like to know what happens in the absence of a clear sep- 
aration of scales between hard and soft modes or when the initial fields have large 
amplitudes which is seemingly the case in present day experiments. Such nonper- 
turbative studies may be based on simulations of lattice gauge theory using real- 
time stochastic quantization techniques @, or classical-statistical methods IfTUII in 
the presence of sufficiently large occupation numbers. Classical-statistical simula- 
tions on a lattice have been extremely successful in the past in describing the non- 
linear scalar inflaton dynamics in the early universe in the presence of parametric 
resonance or spinodal instabilities after inflation 0Tl[T2"in"3l . These instabilities 
lead to exponential growth of occupation numbers in long wavelength modes on 
time scales much shorter than the asymptotic thermal equilibration time. Charac- 
teristic far-from-equilibrium phenomena, such as an early prethermalization of the 
equation of state, have been quantitatively studied in that context B2l Q2|. Com- 
parisons with calculations including quantum corrections demonstrate the accu- 
rate description of nonequilibrium quantum field dynamics by classical-statistical 
methods for not too late times, i.e. before the thermalization characterized by 
Bose-Einstein statistics sets in 0"3l . 

In this paper we compute nonequilibrium dynamics for classical-statistical SU(2) 
pure gauge theory on a lattice. The theory can be computed without further 
approximations using standard techniques by numerical integration and Monte 
Carlo sampling over initial conditions for given normalized initial probability 
functional 031. We consider anisotropic initial conditions with large occupation 
numbers in the transverse plane with the characteristic scale A that may be iden- 
tified with the saturation scale ~ Q s 031. Its value is varied in a range 0.5 < 
A/e 1/4 < 2.5, where e denotes the initial energy density. This is used to inves- 
tigate the very early stages of the thermalization process (isotropization, prether- 
malization) such that a static geometry, i.e. no expansion may be considered. We 
find Weibel or "primary" instabilities with characteristic inverse growth rates sim- 
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ilar to those obtained from previous treatments employing anisotropic distribu- 
tions of hard modes (particles) in the weak coupling limit. We observe "sec- 
ondary" growth rates for higher-momentum modes reaching substantially larger 
values. The secondaries are fluctuation induced and we analyze them in terms 
of resummed loop diagrams beyond the hard-loop approximation. We find that 
a coarse grained pressure isotropizes "bottom-up" with a characteristic inverse 
growth rate of y" 1 - 1 -2fm/c for coarse graining momentum scales of p < 1 GeV 
choosing a value for RHIC of e = 30GeV/fm 3 . We finally analyze the nonequi- 
librium spatial Wilson loop which is found to exhibit an area law. This is used 
to compute the time evolution of a "spatial string tension" which is shown to 
be non-vanishing and isotropic on a similar characteristic time scale as the low- 
momentum coarse-grained pressure. 

The paper is organized as follows. In section [21 we outline our approach and 
the numerical implementation. Subsequently, we analyze the characteristic time 
scales associated to primary and secondary growth rates in section [31 Section 
[4] discusses isotropization of a coarse grained pressure and section [5] computes 
nonequilibrium dynamics of large Wilson loops. We conclude with a discussion in 
section[6]and verify volume and cutoff independence of our results in an appendix. 



2 Classical-statistical gauge field theory on a lattice 

We consider the Wilsonian lattice action for SU(N) gauge theory in Minkowski 
space-time 

S[U] = -^o2Z{2T7l( Tr ^ + Tr ^o,)-l} 

^ZZl^^^' + Trf/g-l}. (1) 

It is described in terms of the plaquette variable 

Ux, MV =U x , M U x+fl , v ul 9M Ul v , (2) 

where U\ Vil = U x ^ v . Here U XifJ is the parallel transporter associated with the link 
from the neighboring lattice point x + p. to the point x in the direction of the lattice 
axis ji = 0, 1, 2, 3 with U X4l = E/jLj . The parameters 

2yTrl 2Trl 
J3 = JL ^~ , J3 S = —r- , (3) 

g Ss7 
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contain the anisotropy parameter y = aja t , where a s denotes the spatial and a t 
the temporal lattice spacings, and we will consider g = g s = g as the coupling 
constant of the lattice theory. 

In the following we specify to SU(2) as the gauge group. The dynamics is solved 
in temporal axial gauge, i. e. U x fi = 1 and the numerical treatment is similar as 
in Ref. ifToll . Varying the action (QQ) w.r.t. the spatial link variables U x j yields the 
leapfrog-type equations of motion 

E b j(t + a t ,x) = E%x) + 1 ^—Y^xi^U^U^U} U\ k ) 
+ Tt(o b U xi ti] . Ul U) . U (xhk )\, 

v X ' J (x+j-k),k (x-k),j (*-*)>*' J ' 

where cf 1 {a = 1, 2, 3) are the Pauli matrices. The electric fields are given by IfTTTl 

E J (X) = - l oVTo ^(U^jUiy) = ~-±- Tr (Urtj**) ■ (5) 
2 a s a t g ,J I a s a x g 

Thus E h Xt + a t , x) can be computed from the electric fields and the spatial link 
variables at time t. Then the link variables at t + a t are given by 

U ix+t)J = | (l - (a s a t g) 2 ^ E)(t + a t , x) 2 ) 1 ' 2 t + ia s a t gE)(t + a t , x)(A U xj . 

(6) 

Varying the action ([T]) w.r.t. to a temporal link gives the Gauss constraint 

3 

J][E b j (x)-Ul_ ]j E b j ( X -j)U x _ lj ]=0. (7) 

;=i 

The coupling g can be scaled out of the equations of motion and we will set 
g = 1 for the simulations. While g is irrelevant for the classical dynamics it 
appears when calculating observables such as pressure as described below. The 
relation between the link variables and Aj is given by Uj(x) = e lgsUsA ^ where 
Aj(x) = A^xjcr" and, accordingly, we use 

A)(t, x) - --i- TT(o*Uj(t, x)) . (8) 

We will often consider the Fourier transform w.r.t. the spatial coordinates, i.e. 
A b j(t,p). 

The initial conditions are given in terms of a normalized initial Gaussian prob- 
ability functional P[A(t = 0), E(t = 0)] from which correlation functions, such as 
the two-point function, are obtained as 

<A?(*)A;O0> = J DA(0)DE(0)P[A(0),E(0)]^(x)A b j(y), (9) 



where the measure indicates integration over classical phase space. The initial 
Gaussian probability functional is chosen such that 

2 2 2 

<l4<' = 0.P)r-> = ^4 3Z; expj-^-^|, (10) 

and we consider E a . = -A" = at t = fulfilling the Gauss constraint ©, 
which ensures that the latter is respected at all times. In practice, we compute 
expectation values by solving the equations of motion © and © for a set of 
initial configurations sampled according to Eq. (PTOl) and then average over the 
results from the individual runs. The margins of errors that we quote are the 
statistical errors. The widths appearing in Eq. (fTUl) are typically chosen such that 
A :» A, and A z is set to a sufficiently small value to have an effectively 5{p z )- 
like distribution on the lattice. We also discuss the effects of broadening the initial 
longitudinal distribution. Here A determines the typical transverse momenta of the 
gluons and may be associated with the saturation scale of order GeV. The constant 
C appearing in (ITOT) is adjusted to obtain a given energy density. For numerical 
reasons, we initialize all Fourier coefficients whose absolute value is less than a 
small lower bound with a little noise term. We have checked that our results are 
insensitive to the amplitude of the noise. 

The energy density in lattice units is determined by the Hamiltonian IfTTll corre- 
sponding to (QQ) which reads 



^ x) s h B E p - Tr u *»a + Z t 2 - Tr u «d 

[ 1 j j<k 



(11) 



For the conversion to physical units we will fix the lattice spacing a s from the 
relation between the physical energy density and its lattice analogue according to 

e = e-a- A . (12) 

If g is taken to be different from one then the lattice spacing is altered by a factor 
1/ V? ' wn i cn follows from Eq. (fTT|) . and the values for a will not be altered 
significantly as long as g ~ 0(1). 

We have checked that the results, which are presented in the following, are 
insensitive to the employed volume and lattice spacings. Details are given in an 
appendix. 



3 Characteristic early time scales 

Instabilities provide the parametrically fastest nonequilibrium processes, which 
may help to understand the apparent success of hydrodynamics in describing 
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Figure 1: (Color online) Fourier coefficients of the squared modulus of the gauge 
field versus time for three different momenta parallel to the z-axis. The low- 
momentum modes exhibit exponential growth from the beginning. This is fol- 
lowed by a secondary stage that sets in later but with significantly larger growth 
rates for modes. The curves correspond to the momenta in the legend from top to 
bottom in the same order. 



RHIC data at times of less than about 1 fm/c. A crucial assumption for the use 
of ideal hydrodynamic models is that the stress tensor, in the local rest frame at 
some point in the system, is nearly diagonal, T t j ~ 6jj, with some equation of state 
relating the pressure to the energy density. This is just a statement of isotropy 
which may be faster achieved than complete local thermal equilibrium [4j. In 
particular, hydrodynamics is seen to provide a description of the data only for not 
too high characteristic momenta of order GeV [QQ|. Therefore, an isotropic coarse 
grained pressure, which does not take into account anisotropics on distance scales 
smaller than some characteristic inverse coarse graining momentum may be suffi- 
cient. Our following results seem to support this possibility. Here we concentrate 
on the very early stages of the dynamics, which is supposed to be quantitative 
only for sufficiently fast processes such that the expansion of the system can be 
neglected. We also emphasize that studying the late-time dynamics is beyond any 
classical-statistical approximation, which cannot explain the approach to thermal 
equilibrium characterized by Bose-Einstein distributions. For sufficiently short 
times and large range of highly occupied momentum modes, however, it is ex- 
pected to provide a quantitative description of the nonperturbative dynamics lfT3ll . 
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Figure 2: Primary and secondary growth rates for \A(t, p)| 2 as a function of p z as 
measured on 96 3 - and 128 3 - lattices using the same lattice spacing. 



3.1 Primary and secondary growth rates 

In Fig. \T\ we plot the color- averaged squared modulus of three different Fourier co- 
efficients of the gauge field modes A b x (t, p) as a function of the difference At = t-t 
to the initial time to, normalized by the corresponding field values at t = to- All 
values are given in appropriate units of the initial energy density e. Here \A(t, p)| 2 
may be associated to particle numbers lfT8ll . The momenta p are chosen parallel 
to the axis of anisotropy, i.e. the z-axis. These modes exhibit the fastest changes 
in qualitative agreement with analytical results in the weak coupling limit for the 
growth of Weibel instabilities in non-Abelian gauge theories [;3J. The plotted low- 
momentum modes clearly show exponential growth starting at the very beginning 
of the simulation. In contrast to these "primary" instabilities operating at low 
momenta, one observes from Fig. \T\ that momentum modes at sufficiently high 
momenta do not grow initially. The higher wave number modes typically exhibit 
exponential growth at a secondary stage that sets in later but with a significantly 
larger growth rate. In particular, for not too high momenta the primary growth 
abruptly changes and a secondary, larger growth rate takes over as exemplified 
for p = 0.8 £ 1/4 in Fig. [TJ The observation of primary and secondary growth 
rates closely resembles findings in the context of early universe preheating dy- 
namics in self-interacting scalar inflaton models following parametric resonance 
or spinodal instabilities lfT9l . Similar to the analysis for inflaton dynamics, we will 
demonstrate below in section lX2l fhat the secondaries arise from fluctuation effects 
induced by the growth in the lower momentum modes by taking into account re- 

'The data is an average of 100 runs on a 64 3 - lattice. 
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Figure 3: Time scales associated with the maximum primary growth rate as a 
function of the transverse width A in units of e ~ 1/4 . 



summed loop diagrams beyond the hard-loop approximation. This growth satu- 
rates when all loop diagrams are of order one corresponding to a highly non-linear 
evolution^ For the calculation for Fig.Q], which employs A/e 1/4 = 1, saturation 
occurs around t ~ 160e -1/4 . 

From Fig. [Hone also concludes that the same phenomena could be observed if 
the simulation started at a later time than to. Choosing the zero point of the time 
axis just changes the details of the initial conditions. For instance, even at times 
At = 50 6 1/4 or 100e~ 1/A the anisotropy is still substantial. As an indicator one 
may consider the bulk anisotropy parameter g(t) defined as 



At ~ 120e~ 1/4 would equally well qualify as initial conditions with a very high de- 
gree of anisotropy. In order to extract characteristic time scales for isotropization 
we, therefore, focus on the inverse growth rates. 

Fig. [2] displays the momentum dependence of the growth rates for \A(t, p)| 2 ob- 
tained from a fit to an exponential. According to the above discussion this is done 
separately for the primary and secondary growth rates. One observes that the sec- 
ondaries occurring at higher momenta exhibit growth rates which are typically 

2 This corresponds to (AA) ~ 1 /g 2 in a description where the coupling would not be scaled out 
as described in Sec. [2] See also the discussion in Ref. lfT4l[l9l . 
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Figure 4: Time evolution of | A(t, p) | 2 from initial conditions with a slightly broad- 
ened distribution in the longitudinal direction with A/A z = 25. The curves corre- 
spond to the momenta in the legend from top to bottom in the same order. 



larger than the highest primary rate by a factor of approximately three to four. 
The dependence of the results on the initial width A in the transverse plane can be 
inferred from Fig. [3j where the inverse of the maximum primary growth rate as a 
function of A/e 1/4 is plotted. 

For the conversion to physical units we will consider two different values for the 
energy density of about 1 GeV/fm 3 and 30 GeV/fm 3 , equivalent to approximately 
(300 MeV) 4 - (1.5/fm) 4 and (700 MeV) 4 =* (3.5/fm) 4 . For RHIC experiments the 
former estimate is expected to be too low, while the latter can be still appropri- 
ate lTT5l . According to Fig. [31 the inverse of the maximum primary growth rate for 
| A(t, p) | 2 is y"^ pr - 4 e~ 1/4 . In physical units this corresponds to 

y m ! x . pr . - l.lfm/c (6 = 30 GeV/fm 3 ), (14) 
y m i xm , - 2.7fm/c (e = 1 GeV/fm 3 ). (15) 

In order to achieve growth by a few orders of magnitude in the low-momentum 
occupation numbers a time of several y m !, x has to elapse. One may then ask 
what energy density would be required to get characteristic inverse growth rates 
of about O.lfm/c and we obtain: 

TrnLp, - 0.1fm/c (6 = 300TeV/fm 3 !) (16) 

The situation changes somewhat if one considers the characteristic time scales 
associated with the secondary growth rates. For these we find that y"^ sec 
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Figure 5: Growth rates of | A(t, p) | 2 as a function of p z corresponding to Fig. HI 



1.3 e 1/4 , which in physical units corresponds to 

TmLsec. - 0.4fm/c (6 = 30GeV/fm 3 ), (17) 
Tmix.sec. - 0.8fm/c (e = 1 GeV/fm 3 ). (18) 

At first sight, these seem to be easier to reconcile with the envisioned fast isotropi- 
zation. However, even though secondaries can reach much higher growth rates 
their appearance can be delayed. According to the results in Fig. \T\ the secon- 
daries only affect the isotropization process at early times if the initial conditions 
happen to be similar to the gauge field configuration at At - 100 - 120e 1/4 be- 
cause otherwise they would set in too late. So far, we have initialized our runs 
using a 6{p z )-like distribution at t = t . Other initial conditions can provide for a 
more pronounced role of secondaries. An example is given in Fig.Hl which shows 
the time evolution of | A(t, p) | 2 starting from the same Gaussian initial condition 
as given in Eq. (flOl) . however, for a slightly broadened distribution in the longi- 
tudinal direction with A/A z = 25. Indeed, one sees that from the very beginning 
growth rates comparable to the secondary instabilities observed in Fig.[T]occur. In 
Fig.[5]these rates are given as a function of p z . Comparison with Fig.[2]for given 
momenta shows that the growth is even faster by a factor of approximately two. 
This would lead to somewhat shorter time scales than presented in Eqs. (fT4l) -(PT8l. 
The analysis of (primary) growth rates using hard loop approximations leads to the 
expectation that the rates decrease with decreasing A/A z , i.e. the less anisotropic 
the initial conditions are. The above results show that (secondary) growth rates 
can, however, increase with decreasing A/A-. Their understanding requires to go 
beyond hard loop approximations, which is done in the following. 
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(a) 



(b) 



Figure 6: Loop diagrams (a) UfZ (t, p) given in Eq. ([2"2l and (b) II ^ (t, p) given 
in Eq. (|23T) . The propagator lines correspond to the exact classical-statistical prop- 
agator, while the vertices are the perturbative ones. 



3.2 The impact of fluctuations: Diagrammatic analysis 

We have seen above that low momentum modes of | A(t, p) | 2 show an exponen- 
tial increase starting at the very beginning of the simulation, while for sufficiently 
large initial anisotropy higher momentum modes exhibit exponential growth at a 
secondary stage with rates several times larger than primary ones. Following the 
discussion of qualitatively similar findings in Refs. | fT4l [T9ll for self-interacting 
scalar preheating dynamics, we analyze the appearance of secondaries in terms 
of resummed self-energy diagrams. These diagrams appear as source terms in 
the Schwinger-Dyson equations for propagators. They naturally arise using two- 
particle irreducible (2PI) effective action techniques for the description of nonequi- 
librium dynamics of correlators such as the equal-time two-point function |fl4| 



Since we work in temporal axial gauge, i.e. A a Q {x) = 0, it follows that F^ = = 

pab A 

Low-order diagrams in a loop expansion of the self-energy contributing to the 
correlator (fT9l) are displayed in Fig.[6l We emphasize that the diagrams are com- 
puted using the exact classical-statistical propagator obtained from simulations 
and perturbative vertices 031 |20l . The perturbative three- and four-gluon vertices 
are given as 



F%(t,x-j) = {Afi,x)A$(t,j)). 



(19) 




(k, p, q) = ge abc (g^(k - pf + £"(p - qf + ^(q - k) v } 



(20) 
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and 



KZ d = -i8 2 {e ah fe cdf (g^g v<T ~ ^V P ) 

+ eacfeufbTfr ~ S"V P ) + e adf e bcf (g^r ~ ffg™)} ■ 

In the above vertices the Lorentz indices are restricted to run over spatial directions 
only (i. e. // = {1,2,3}) since temporal components only appear together with 
vanishing and F"^. The corresponding contribution of the one loop diagram 
Fig.[6ta) reads 

K/ 1 \t, P) = \ J H { VZ C (~V - q, P, q)V$.(p + q, "P, -q) ^ 

x F>,p + q)F^,-q)(. 

The two-loop diagram Fig.HJb) is evaluated in coordinate space according to 

rC (2 V, x - y) = i ValV^JCf, x - y)F/f (f, x - y)F£fc x - y) (23) 

and then Fourier transformed. 

Fig. |7] shows the results of such a calculation where all external color indices 
have been set to one and Lorentz indices were chosen in the transverse plane, 
i.e. perpendicular to the axis of anisotropy. The upper graph displays F|}(f, p) 
for a given momentum as a function of time, similar to Fig. [H The momentum 
p z /e 1/4 = 1 is chosen such that the mode exhibits a primary growth rate initially 
which abruptly changes at a later time of about Ate 1/4 = 100 to a secondary rate. 
The lower graph of Fig. [7] gives the absolute values of the loop diagrams Eqs. (1231) 
and (1221) normalized to the correlation function Eq. (fT9l) . The diagrams are evalu- 
ated for the same momentum as the gauge fields in the upper graph. One observes 
that initially these loop corrections are small compared to the size of the correla- 
tor (TT9l) in units of e. Accordingly, their contribution to the time evolution of the 
correlator may be neglected and primary growth rates appear as a consequence 
of Weibel instabilities as discussed in great detail in the literature 0. However, 
due to these instabilities loop corrections grow exponentially and faster than the 
correlator itself. The faster growth is seen to stop at about At = 100 e~ 1/4 . This is 
precisely the time where the secondary growth rate for the correlator sets in. The 
faster growth of the loop corrections originates from the fact that the exponentially 
growing correlation modes F(t, p) appear multiple times as propagators in the loop 
diagrams such as shown in Fig. [6] For the case of the scalar preheating dynam- 
ics analyzed in Refs. IfFfl it was shown that secondary growth rates indeed 
equal multiples of primary ones. This was possible because primary growth was 
sharply peaked around certain momenta in the scalar models at weak coupling, 
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Figure 7: The upper graph displays the correlation function F(t, p) defined in 
Eq. (fl9l) as a function of time for given momentum p z /e i/4 = 1. The lower graph 
shows for the corresponding momentum the relative contributions of the diagrams 
p)/i 7 (f, p) and n (2) (Y, p)/F(t, p) to the correlation function. The upper curve 
gives the one-loop contribution, the lower one is the setting-sun contribution. 



and the loop integrals could be estimated in a saddle point approximation around 
the dominating loop-momentum. For the SU(2) gauge theory considered here pri- 
maries do not appear in a sharply peaked momentum range as can be observed 
from Fig. 121 Accordingly, we do not find that secondary growth rates correspond 
to simple multiples of primary ones. 

In section [3TTI we have seen that for large anisotropics, corresponding to S(p z )- 
like initial distributions (PTOl) . sizeable secondaries only appear at later times and 
not from the beginning of the simulation. This can be understood from the fact 
that for these initial conditions the loop corrections as in Fig. [6] are negligible at 
initial time for longitudinal momenta. In contrast, broadening of the initial dis- 
tribution in the longitudinal direction can lead to sizeable initial loop corrections. 
Accordingly, we find for such initial conditions larger growth rates already at ear- 
liest times as demonstrated in Figs. [4] and [51 However, from comparing Fig. [4] 
with Fig. [T]one observes that the exponential behavior is less developed the less 
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Figure 8: Time evolution of longitudinal "pressure", i.e. the spatially Fourier 
transformed longitudinal component of the energy-momentum tensor T 33 (f, p) 
normalized to initial values for various momenta p z . 



anisotropic the initial conditions. As a consequence, we find that characteristic re- 
sults for the isotropization of observables such as pressure are rather insensitive to 
the details of the considered initial conditions which is discussed in the following 
section. 



4 Bottom-up isotropization 

We compute the time evolution of the spatial components of the energy-momentum 
tensor with initial conditions as employed for Figs.[T]-0 Its diagonal components 
yield the pressure in thermal equilibrium and we will refer to Tn and T22 as the 
time-dependent transverse pressure. Accordingly, r 33 will be called the longitu- 
dinal pressure. The diagonal components of T, 7 can be computed according to 

Txix) = - ^ (y 2 ((l - l - Tr U m ) + (1 - i Tr U xfi2 ) - (1 - \ Tr U xfi3 )) 

+ (1 - X - Tr t/,, 2 3 ) + (1 - X - Tr U xM ) - (1 - X - Tr U x , 12 )} , 

(24) 

and Tn and T 22 can be obtained by permuting the indices respectively. We have 
verified that non-diagonal components of the energy-momentum tensor are neg- 
ligible at all times for the class of initial conditions considered here. We then 
Fourier transform w.r.t. the spatial components and consider the momentum modes 
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Figure 9: Snapshots of the spectrum of r 33 (7, p)/T 33 (Y = 0, p) at various times. 



of the longitudinal and transverse pressure as a function of time. Fig. [H] shows the 
longitudinal pressure modes for four different momenta in the z-direction, normal- 
ized to the respective initial values. One observes that the longitudinal pressure - 
which is a gauge invariant quantity - exhibits instabilities in a similar way as does 
the gauge field. The same qualitative picture regarding time scales, primary and 
secondary instabilities emerges from Fig.[H]as from Fig.[TJ Quantitatively, a closer 
look reveals that the growth rates of modes of T 33 are smaller by a factor of about 
two than the corresponding modes in | A(p) | 2 with the same wave number. This 
can be understood from the fact that the longitudinal pressure contains contribu- 
tions 2 q Mt % p)A(Y, q - p), which are dominated before saturation by momenta 
where q - p lies in the transverse plane. Since the transverse modes do not exhibit 
instabilities, there is only one exponentially growing factor appearing in each of 
the dominant contributions. Fig. |9] shows the same quantity as in Fig. [8l however, 
as a function of longitudinal momentum. Plotted are snapshots of the spectrum 
at fixed times versus the absolute value of the momentum. One observes a band 
of low-momentum modes that are unstable at initial times. At the time when the 
secondary instabilities set in also a broad regime of higher modes suddenly starts 
to grow, similar to what has been described in Ref. [EH and in the second paper 
of Ref. 0. 

The momentum scale of the Fourier transformed spatial components of the 
energy-momentum tensor may be viewed as a characteristic scale of a coarse 
grained pressure. One can then ask the question about a characteristic coarse 
graining scale up to which fast isotropization may occur. Fig. [lOl shows the ratio 
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Figure 10: The ratio of modes of longitudinal pressure and transversal pressure 
(1251) as a function of time for different momenta. One observes that pressure 
isotropizes "bottom-up" and complete isotropization does not occur on short time 
scales. 



of longitudinal and transverse modes of pressure, i.e 

Tyi(t, Ph) 



(25) 

I ph 1=1 p± \=p 



for three different momenta as a function of time. One observes that for low- 
momentum modes the ratio turns to one rather quickly. High momentum modes 
only become isotropic at far later times. Complete isotropization does not occur 
on short time scales. The system becomes only effectively isotropic when consid- 
ered on sufficiently large characteristic length scales. This may be used to define 
an isotropic coarse grained pressure by averaging over a volume with a charac- 
teristic length scale given by the inverse momentum. For instance, requiring that 
Pl(p)/Pt(p) £ 0.6 at the time when the instabilities saturate (At = 160e~ 1/4 in 
Fig. [[]) we find isotropization up to momenta p z < 1.4e 1/4 , which corresponds to 

p z < lGeV (e = 30 GeV/fm 3 ) , (26) 

p z < 400MeV (6 = 1 GeV/fm 3 ). (27) 

Hydrodynamic simulations indeed assume locally isotropic pressure only on suf- 
ficiently large characteristic length scales [1J. However, even for e = 30 GeV/fm 3 
the maximum coarse graining momentum we find is somewhat small. These re- 
sults are obtained for the (5(p.)-like initial conditions as explained above, and in 
view of the results of section [3] one might expect that broadening the initial dis- 
tribution in the longitudinal direction improves the situation and leads to earlier 
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isotropization of a broader momentum range. However, we find that the maximum 
mode that has a ratio of Pl(~P)/Pt(v) 0.6 at saturation has a momentum of about 
p z - 1.5e 1/4 for the less anisotropic initial conditions employed for Figs. |4] and [51 
Despite the fact that growth rates are twice as large in this case, the exponential 
behavior is less developed the less anisotropic the initial conditions. We conclude 
that the estimates (1271) are rather insensitive to the details of the considered initial 
conditions. 

5 Wilson loops 

Wilson loops as traces of parallel transporters around rectangular loops provide 
an important means to describe confinement properties of non-Abelian gauge the- 
ories. In vacuum large Wilson loops with side length R and R obey the area law 

W(R,R) ~ e~ KRR (28) 

for R,R — » oo. Here k denotes the string tension and a non- vanishing value can 
be related to static quark confinement for loops including a time separation on a 
Euclidean space-time lattice. For the latter it is known that spatial Wilson loops, 
where R and R both denote spatial distances, give the same value for k as the 
temporal ones and that the "spatial string tension" depends rather weakly on the 
temperature 11221 . 

On a real-time lattice the nonequilibrium time evolution of spatial Wilson loops 
can be measured and the corresponding spatial string tension can be extracted 
according to Eq. (|28T) . We measure the time evolution of large Wilson loops in 
the transverse plane, i.e. with R and R measuring distances in the xy-plane, and 
loops where one side length measures distance in the z-direction. We will call the 
former transverse Wilson loops, W T , and the latter longitudinal Wilson loops, W L , 
accordingly. We emphasize that since we are considering classical-statistical time 
evolution we do not take into account quantum fluctuations. This can only be a 
good approximation if classical-statistical fluctuations dominate. While this is less 
problematic for bulk quantities such as pressure, which is not dominated by low- 
momentum fluctuations, this is in general not justified for a quantity such as the 
string tension that is sensitive to long-ranged fluctuations related to confinement. 
Accordingly, the longitudinal Wilson loops do not show an area law at initial time 
for the class of initial conditions considered here because of the low occupation 
numbers in the longitudinal direction. 

It is remarkable, however, that the transverse Wilson loops exhibit an area law as 
a consequence of large classical-statistical fluctuations in the presence of the high 
initial transverse occupation numbers. This is demonstrated in Fig. [CD where the 
values for Wl and Wj are given at various times as a function of the area R x R 
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Figure 1 1 : Values of spatial Wilson loops in the transverse plane (W T ) and Wilson 
loops including the longitudinal direction (W L ) as a function of the loop area in 
lattice units. Shown are snapshots at increasing time going to the right. Initially, 
Wl exhibits no area law, while the transverse Wilson loop is characterized by 
a non-vanishing spatial string tension k as explained in the text. One observes 
W L /W T - 1 at the time of saturation of the exponential growth of gauge field 
fluctuations and the anisotropy is removed. 



in lattice units. We evaluate a variety of rectangular loops with one side length 
R = 2, 4 or 6, while the other side length is varied in the range R e {2, 4, 20} 
and square-shaped Wilson loops with R = R. One observes that initially W L ex- 
hibits no area law. In contrast, the values of the transverse Wilson loops follow an 
area law to very good accuracy. The initial conditions are as described in section[2] 
for a £(p z )-like initial distribution and we employ a lattice energy density e = 0.05 
according to Eq. (fT2l) and a transverse width A = 0.5 in lattice units on a 64 3 
lattice. The figures to the right show the same quantities at later times. The lon- 
gitudinal Wilson loop is seen to develop an area law and at the time t 160e^ 1/4 , 
after the exponential growth of gauge field fluctuations saturates, all Wilson loops 
follow the same area law and the anisotropy is removed with Wl/Wt — 1. We 
conclude that the isotropization of the spatial Wilson loops occurs on a similar 
characteristic time scale as the low-momentum coarse-grained pressure. 

Though the string tension extracted from the transversal Wilson loops is non- 
vanishing from the beginning, it changes with time. In FigH2]the corresponding 
value of /c 1/2 is plotted as a function of time in units of the energy density e. For 
the considered initial conditions one observes an initial decrease with time, which 
tends to a non-zero value /c 1/2 /e 1/4 ^ 0.14 at late times. We verify the expectation 
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Figure 12: The time dependence of the spatial string tension from Wilson loops 
in the transverse plane. 



that quantities such as pressure are not sensitive to this scale. As described in 
sections [3] and [5] the pressure as well as the growth rates only depend on A/e 1/4 for 
a given time or momentum. 

The result of a non-vanishing scale as defined by Eq. (1281) for spatial Wilson 
loops may seem surprising in a classical-statistical theory at first sight. How- 
ever, our nonequilibrium results relate well to known behavior in thermal equilib- 
rium at high temperature. The low-momentum behavior of a 3 + 1 -dimensional 
bosonic quantum theory at sufficiently high-temperature is described in terms of a 
3-dimensional classical-statistical theory by virtue of dimensional reduction, and 
the 3-dimensional non-Abelian gauge theory is known to exhibit an area law in 
thermal equilibrium [231. O ur results may be viewed as nonequilibrium general- 
ization of this. 



6 Discussion 

Real-time lattice QCD is a crucial tool for the quantitative understanding of the 
nonperturbative thermalization processes in heavy-ion collisions. Late-time dy- 
namics requires calculation of nonequilibrium dynamics in the full quantum the- 
ory, which is still a formidable task. Important aspects of the short-time dynamics 
in the presence of high gluon occupation numbers, however, may be calculated 
using the classical-statistical field theory limit without further approximations. In 
particular, for a large class of initial conditions there is no need to model high 
momentum degrees of freedom as particles or to restrict the discussion to the very 
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weak coupling limit. In this work we have considered classical-statistical SU(2) 
gauge theory for anisotropic initial conditions with high occupation numbers in 
the transverse plane on a characteristic scale ~ A. We observe Weibel or primary 
instabilities with growth rates similar to those obtained from previous treatments 
using a variety of setups and approximations (31 HI [51 [6). We observe secondary 
growth rates for higher-momentum modes reaching substantially larger values. 
The secondaries are fluctuation induced and reflect strongly non-linear dynamics. 

Two types of gauge invariant observables have been calculated: Firstly the 
time evolution of the spatial part of the energy momentum tensor, which is re- 
lated to pressure. Its Fourier transform may be considered as a scale dependent 
pressure with a characteristic coarse-graining momentum scale. We observe that 
plasma instabilities lead to a rapid "bottom-up" isotropization up to momenta of 
about lGeV for an energy energy density of e = 30GeV/fm 3 . The characteris- 
tic time scale is given by the inverse growth rate of about 1 - 2fm/c. Complete 
isotropization does not occur on short time scales. Therefore it is important to 
note that an effectively isotropic pressure for a low-momentum range of character- 
istic momenta is a crucial ingredient for the use of hydrodynamics, and complete 
isotropization may not be necessary to account for its successful applications to 
describe RHIC data. 

The pressure is a quantity which is not sensitive to long-ranged fluctuations 
in contrast to the second type of observables we compute: The nonequilibrium 
spatial Wilson loop is shown to become isotropic on similar characteristic time 
scales as the coarse grained pressure for low momenta. At this time all large 
spatial Wilson loops are found to obey an area law, which is used to compute a 
non-vanishing scale with the dimension of a string tension. The generation of 
a non-vanishing scale demonstrates the non-perturbative nature of the classical- 
statistical limit of the corresponding quantum theory. We have verified that the 
pressure and growth rates are insensitive to this scale, which is required for weak- 
coupling descriptions of these quantities. 

The extension to SU(3) pure gauge theory is straightforward, though compu- 
tationally more expensive, in order to verify that the changes are indeed minor 
as expected. Since fermions cannot have large occupation numbers, their contri- 
butions to the low-momentum dynamics is expected to be small and the results 
should give a quantitative account of the short-time physics. 

We thank Adrian Dumitru, Carsten Greiner, Mike Strickland and Zhe Xu for 
very interesting discussions, and Mike for re-calculating some of the results pre- 
sented here using his setup. This work is supported in part by the BMBF grant 
06DA267, and by the DFG under contract SFB634. 
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Figure 13: (Color online) Fourier coefficients of the gauge fields as in Fig. Q] 
computed for three different lattice volumes. The left panel shows the momentum 
p = 0.5 e 1/4 while the right panel is for p = 0.8 e 1/4 . While the curves computed 
on 64 3 - and 128 3 - lattices coincide well, the one from the 32 3 - lattice (red) lies a 
bit lower. 



A Appendix 

In this appendix we demonstrate the insensitivity of our results under variations of 
volume or spatial lattice distance. All simulations were done with g = 1 and the 
initial configurations were generated from the momentum space probability distri- 
bution Eq. (flOl) using the same values for the widths and fixing the normalization 
such that the energy densities agree. 

First, we check for possible volume dependences. In Fig. [T3]fhe time evolution 
of the squared modulus of the gauge field Fourier modes is plotted similarly as 
in Fig. [Has obtained from lattices of 32 3 , 64 3 and 128 3 sites. On all lattices, we 
have chosen modes parallel to the z- axis with equal magnitude p 0.5 6 1/4 and 
p - 0.8 6 1/4 for the left and the right graph of Fig. \\3\ respectively. Volume 
dependences can hardly be observed for lattices of size 64 3 or larger and also on 
the smallest lattice (32 3 ) only minor deviations occur. 

Since the UV part of the spectrum is not excited in the beginning, the depen- 
dence on the UV cutoff or, equivalently, the lattice spacing should be even milder 
than the volume dependence. In Fig. [14] the time evolution of Fourier modes is 
plotted as obtained from simulations on lattices of different lattice spacings but 
same volume. We have chosen the first and the second non-zero momentum of 
the 32 3 - lattice in the z- direction. The lattice spacings are a = 0.25,0.5, 1,2 in 
units of the lattice spacing used for the results presented in the main body of this 
paper. One observes from Fig. [14] that the results are insensitive to UV cutoff 
changes if the lattice spacing employed in the rest of this paper or smaller is used. 
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Figure 14: (Color online) The time evolution of the first and second non-zero 
mode parallel to z-axis plotted for systems with different lattice spacings. Time is 
measured in units of the lattice spacing of the 32 3 - lattice. The left graph employs 
the momentum p = 0.20 /a^ while the right graph uses p = 0.39/^32. The order- 
ing of the legend corresponds to the ordering of the curves from top to bottom at 
t = 0. 



22 



References 

[1] U. W. Heinz, AIP Conf. Proc. 739 (2005) 163. 

[2] J. Berges, S. Borsanyi and C. Wetterich, Phys. Rev. Lett. 93 (2004) 142002. 

[3] S. Mrowczyriski, Phys. Lett. B 214 (1988) 587; P. Romatschke and 
M. Strickland, Phys. Rev. D 68 (2003) 036004; P. Arnold, J. Lenaghan and 
G. D. Moore, JHEP 08 (2003) 002; P. Romatschke and M. Strickland, Phys. 
Rev. D 68 (2003) 036004; S. Mrowczyriski, A. Rebhan and M. Strickland, 
Phys. Rev. D 70 (2004) 025004; A. Rebhan, P. Romatschke and M. Strick- 
land, Phys. Rev. Lett. 94 (2005) 102303; P. Arnold, G. D. Moore and 
L. G. Yaffe, Phys. Rev. D 72 (2005) 054003; A. Rebhan, P. Romatschke 
and M. Strickland, JHEP 09 (2005) 041; P. Romatschke and A. Rebhan, 
Phys. Rev. Lett. 97 (2006) 252301; P. Arnold, G. D. Moore, Phys. Rev. D 73 
(2006) 025006; ibid. 025013; B. Schenke and M. Strickland, Phys. Rev. D 
74 (2006) 065004; B. Schenke, M. Strickland, C. Greiner and M. H. Thoma, 
Phys. Rev. D 73 (2006) 125004. 

[4] P. Arnold, J. Lenaghan, G. D. Moore and L. G. Yaffe, Phys. Rev. Lett. 94 
(2005) 072302. 

[5] A. Dumitru and Y. Nara, Phys. Lett. B 621 (2005) 89; A. Dumitru, Y. Nara 
and M. Strickland, Phys. Rev. D 75 (2007) 025016; D. Bodeker and K. Rum- 
mukainen, JHEP 07 (2007) 022. 

[6] P. Romatschke and R. Venugopalan, Phys. Rev. Lett. 96 (2006) 062302. 
ibid., Phys. Rev. D 74 (2006) 04501 1. 

[7] R. Baier, A. H. Mueller, D. Schiff and D. S. Son, Phys. Lett. B 502 (2001) 
51. 

[8] Z. Xu and C. Greiner, Phys. Rev. C 76 (2007) 02491 1. 

[9] J. Berges and D. Sextv. larXrv:0708.0779l fhep-latl: J. Berges, S. Borsanyi, 
D. Sexty and I. O. Stamatescu, Phys. Rev. D 75 (2007) 045007; J. Berges 
and I. O. Stamatescu, Phys. Rev. Lett. 95 (2005) 202003. 

[10] G. Aarts, G. F. Bonini and C. Wetterich, Phys. Rev. D 63 (2001) 025012. 

[11] S. Yu. Khlebnikov and I.I. Tkachev, Phys. Rev. Lett. 77 (1996) 219; 
T. Prokopec and T.G. Roos, Phys. Rev. D 55 (1997) 3768; I. Tkachev, 
S. Khlebnikov and L. Kofman, A. D. Linde, Phys. Lett. B 440 (1998) 262, 
S. Borsanyi, A. Patkos and D. Sexty, Phys. Rev. D 68 (2003) 063512. 



23 



[12] D. I. Podolsky, G. N. Felder, L. Kofman and M. Peloso, Phys. Rev. D 73 
(2006) 023501; J. F. Dufaux, G. N. Felder, L. Kofman, M. Peloso and 
D. Podolsky, JCAP 0607 (2006) 006. 

[13] G. Aarts and J. Berges, Phys. Rev. Lett. 88 (2002) 041603; A. Arrizabalaga, 
J. Smit and A. Tranberg, JHEP 0410 (2004) 017. J. Berges and T. Gasenzer, 
Phys. Rev. A 76 (2007) 033604. 

[14] J. Berges, AIP Conf. Proc. 739 (2005) 3; arXiv:hep-ph/0409233. 

[15] M. Gyulassy and L. McLerran, Nucl. Phys. A 750 (2005) 30. 

[16] J. Ambj0rn, T. Askgaard, H. Porter, and M. E. Shaposhnikov, Nucl. Phys. B 
353 (1991) 346. 

[17] J. B. Kogut and L. Susskind, Phys. Rev. D 11 (1975) 395. 

[18] J.I. Skullerud, J. Smit and A. Tranberg, JHEP 0308 (2003) 045. 

[19] J. Berges and J. Serreau, Phys. Rev. Lett. 91 (2003) 1 1 1601. 

[20] J. Berges, Phys. Rev. D 70 (2004) 105010. 

[21] P. Arnold and G. D. Moore, Phys. Rev. D 73 (2006) 025006. 

[22] E. Manousakis and J. Polonyi. Phys. Rev. Lett. 58 (1987) 847. 

[23] I. Montvay and G. Miinster, Quantum Fields on a Lattice (Cambridge Uni- 
versity Press, Cambridge, England, 1997). 



24 



